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What is the relation between spin squeezing and entanglement? To clarify this, we derive the 
full set of generalized spin squeezing inequalities for the detection of entanglement. These are 
inequalities for the mean values and variances of the collective angular momentum components Jk- 
They can be used for the experimental detection of entanglement in a system of spin-i particles 
in which the spins cannot be individually addressed. We present various sets of inequalities that 
can detect all entangled states that can be detected based on the knowledge of: (i) the mean values 
and variances of Jk in three orthogonal directions, or (ii) the variances of Jk in three orthogonal 
directions, or (iii) the mean values of j\ in three orthogonal directions or (iv) the mean values 
and variances of Jk in arbitrary directions. We compare our inequalities to known spin squeezing 
entanglement criteria and discuss to which extent spin squeezing is related to entanglement in the 
reduced two-qubit states. Finally, we apply our criteria for the detection of entanglement in spin 
models, showing that they can be used to detect bound entanglement in these systems. 

PACS numbers: 03.65.Ud, 03.67.Mn, 05.50.+q, 42.50.Dv 



I. INTRODUCTION 

Entanglement lies at the heart of many problems in 
quantum mechanics and has attracted an increasing at- 
tention in recent years [l], Entanglement is needed 
in several quantum information processing tasks such as 
teleportation and certain quantum cryptographic proto- 
cols. It also plays an important role in quantum com- 
puting making it possible that quantum computers can 
outperform their classical counterparts for several prob- 
lems such as prime factoring or searching. Moreover, en- 
tangled states and the creation of quantum entanglement 
naturally arise as goals in nowadays quantum control ex- 
periments when studying the non-classical phenomena in 
quantum mechanics. 

When in an experiment entanglement is created, it is 
important to detect it. Thus, in many quantum physics 
experiments the creation of an entangled state is followed 
by measurements. Based on the results of these mea- 
surements, the experimenters conclude that the produced 
state was entangled. However, in many-particle experi- 
ments the possibilities for quantum control are very lim- 
ited. In particular, the particles cannot be individually 
addressed. In such systems, the entanglement can be 
created and detected with collective operations. 
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Spin squeezing is one of the most successful approaches 
for creating quantum entanglement in such systems 

[ii&SaSiEEtilQIQEKp- Ref. i 

defined spin squeezing in analogy with squeezing in quan- 
tum optics: Let us consider an ensemble of N spin- 1 par- 
ticles, and define the observables for the collective angu- 
lar momentum as 

fc=i 
( k) 

for I = x,y,z and where a\ are Pauli matrices. Then, 
the variances of the angular momentum components are 
bounded by the following uncertainty relation 

(AJ z ) 2 (AJ,) 2 >i|(J K )| 2 . (2) 

If (AJ Z ) 2 := (J 2 ) — (J z ) 2 is smaller than the stan- 
dard quantum limit ^ | { } | then the state is called spin 
squeezed (l8| . In practice this means that the mean an- 
gular momentum of the state is large, and in a direc- 
tion orthogonal to it the angular momentum variance is 
small. An alternative and slightly different definition of 
spin squeezing considered the usefulness of spin squeezed 
states for reducing spectroscopic noise or to improve the 
accuracy of atomic clocks 0, [l9| . 

It has already been noted in Ref. 0] that the occur- 
rence of spin squeezing is connected to the correlations 
between the spins. In fact, as shown in Ref. [8j, there 
is an entanglement criterion for the detection of the en- 
tanglement of spin squeezed states: If an iV-qubit state 



2 



violates the inequality 

( AJ *) 2 > 1 (3) 

(Jx) 2 + (Jy) 2 ~ [> 

then the state is entangled (not separable), that is, it can 
not be written as [20] 

k 

where the pk form a probability distribution. 

After this first entanglement criterion several general- 
ized spin squeezing criteria for the detection of entan- 
glement appeared in the literature [U I22L I2II and have 
been used experimentally (24, 25]. In Ref. |22|. a general- 
ized spin squeezing inequality was presented that detects 
entanglement close to many-body spin singlets, such as 
for example, the ground state of an anti-ferromagnetic 
Heisenberg chain. In Refs. [2l| . a generalized spin squeez- 
ing criterion was presented detecting the presence of two- 
qubit entanglement. For symmetric systems, these crite- 
ria are necessary and sufficient. In Ref. (23[, other criteria 
can be found that detect entanglement close to symmet- 
ric Dicke states. All these entanglement conditions were 
obtained using very different approaches. Therefore, one 
may ask: Is there a systematic way of finding all such 
inequalities? Clearly, finding such optimal entanglement 
conditions is a hard task since one can expect that they 
contain complicated nonlincarities. 

In Ref. [2||, we have presented a set of spin squeez- 
ing inequalities for the detection of entanglement. We 
showed that these inequalities are complete, in the sense 
that they can detect all entangled states that can be de- 
tected by the knowledge of (Ji) and (AJ;) 2 for three or- 
thogonal directions I — x,y, z. This completeness means 
the following: A state that is not detected by the inequal- 
ities, cannot be distinguished from a separable state by 
knowing (J;) and (A J;) 2 only. 

In this paper we present extensions of this approach in 
several directions. In Sec. II, we first present a detailed 
derivation of the optimal spin squeezing inequalities from 
Ref. [26[ . Then, we consider the case when only the vari- 
ances (AJ/) 2 (and not the mean values (J;)) are known, 
or when only the mean values (J 2 ) are known. We derive 
the optimal spin squeezing inequalities also for this case. 
In Sec. Ill, we consider the case when (Ji) and (AJ;) 2 
are known not only in three orthogonal directions, but 
in arbitrary directions. In this case, we can reformulate 
the spin squeezing inequalities as inequalities for corre- 
lation and covariance matrices. In Sec. IV, we compare 
our optimal spin squeezing inequalities to other known 
entanglement criteria. In Sec. V, we discuss the issue of 
detecting entanglement of the multi-qubit quantum state 
vs. detecting entanglement in the reduced two-qubit den- 
sity matrix. Finally, in Sec. VI we apply our inequali- 
ties to the investigation of spin models. We have shown 
already in Ref. [26j that the spin squeezing inequalities 




FIG. 1: (a) The polytope of separable states corresponding 
to Eqs. (0 for N = 10 and for J = 0. The origin of the coor- 
dinate system corresponds to a many-body singlet state, (b) 
The same polytope for J = (0,0,4). Note that this polytope 
is a subset of the polytope in (a). 



can detect bound entanglement (a weak form of entangle- 
ment, which is at the heart of many fundamental prob- 
lems in entanglement theory) in such models. Here, we 
present more examples for the applicability of the spin 
squeezing inequalities. 



II. THE OPTIMAL SPIN SQUEEZING 
INEQUALITIES 

Our aim is to characterize the separable states in terms 
of the values of (Ji) and (AJ;) 2 . Note that the knowledge 
of (Ji) and (A J;) 2 is equivalent to the knowledge of (J/) 
and (J 2 )- We now present our main result from Ref. 26]: 
Observation 1. Let us assume that for a physical system 
the values of 

J:=((J x ),(J y },(J z }) (5) 
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and 



K:={{Jl),(Jt),{Jl)) 



(6) 



are known. For separable states all the following inequal- 
ities are fulfilled: 



(4) + (J'y) + <J2) < 



(AJ X ) 2 + (AJ V ) 2 + (AJ 2 



f > f , 

N 



(7a) 
(7b) 



(4) + (Jf)-f < (N-l)(AJ m )\ (7c) 

(iV-l)[(AJ fc ) 2 + (AJ0 2 ] > (J 2 m } + ^ 1 , (7d) 

where k, I, m take all the possible permutations of X, y, z. 
While Eq. J7a[ ) is valid for all quantum states, violation 
of any of Eqs. f 7b|| 7rf| ) implies entanglement. 

Proof. The variance, defined as (AA) 2 := (A 2 ) - {A) 2 , 
is concave in the state, that is, if p = pp\ + (1 —p)p2, then 
{AA) 2 p > p{AA) 2 pi + (1 -p)(AA) 2 2 . Thus, it suffices to 
prove that the inequalities of Observation 1 are satisfied 
by pure product states. Based on the theory of angular 
momentum, inequality Eq. (|7a|) is valid for all quantum 
states and the equality holds for states of the symmetric 
subspace. However, for separable states it can be proved 
easily without this knowledge using that for such states 

da 

(4 4) ^ ) ) + (^ ) ^' ) ) + (^ ) ^ ) )<l- (8) 
For Eq. (|7b[) one first needs that for product states 

(AJ fe ) 2 = f -iE^) 2 (9) 

i 

holds. Then, for a product state one has 

(Aj x ) 2 +(AJ y ) 2 +(Aj z ) 2 = ™-iJ24+yl+4. (io) 




< J > 10 
y 



FIG. 2: The polytope of separable states corresponding to 
Eqs. (T7J) for N = 10 and for J — 0. The points corresponding 
to random separable states fill the polytope. 



and the normalization of the Bloch vector, it follows that 



2)>^£(1 



z 2 ) > 0. 



(14) 



Eq. (|7dp can be proved in a similar way. We have to 
show that 



(N-l) [(AJ k ) 2 + (AJ ; ) 2 ] - (J 2 m ) 



4 — u- 

(15) 

This can be proved by rewriting 3 with the individual 
spin coordinates and using Eq. (|13|) : 



3 = (N-l)[ 



4 E X i 



IE* 



> 



^E 



(1 - x 2 - yj - z 2 ) > 0. 



(16) 



Here .r. 



(ay), and Zi 



crz). Knowing 



that xf + yf + zf < 1, the right hand side of Eq. (jT0|) is 
bounded from below by 

Concerning Eq. (fTcjl . we have to show that 

1){AJ X ) 2 + §- {J 2 ) - (J 2 ) > 0. (11) 

This can be written as 



2) = (iV-l)[f -±$> 2 ] 



= (^-i)[f -lE^i 



ii(E^) 2 + (E' 



Using 



(E^) 2 <^E^ 



(12) 



(13) 



□ 

For any value of J the eight inequalities Eqs. define 
a polytope in the three-dimensional ((J 2 ), (J 2 ), ( J 2 ))- 
space. Observation 1 states that separable states lie 
inside this polytope. The polytope is depicted in 
Figs. [lja,b) for different values for J. It is completely 
characterized by its extremal points. Direct calculation 
shows that the coordinates of the extreme points in the 
((Jx), (Jy), (</ 2 ))-space are 



K. := 



B x :— 



^ - K ((J y ) 2 + {J z )% J + <Jy)\ J + «(^> 2 



(Jy) 2 + (J Z ) 2 N 



N 



N 



- + K (J y f- + K {J z } 



where k := (N — 1)/N. The points A y / Z and B y / Z can be 
obtained in an analogous way. Note that the coordinates 
of the points A k and B k depend nonlinearly on (Jfe). 

One might ask whether all points inside the polytope 
correspond to separable states. This would imply that 
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the criteria of Observation 1 are complete, that is, if the 
inequalities are satisfied, then the first and second mo- 
ments of Jk do not suffice to prove entanglement. In 
other words, it is not possible to find criteria detecting 
more entangled states based on these moments. Due to 
the convexity of the set of separable states, it is enough 
to investigate the extremal points: 

Observation 2. (i) For any value of J there are sepa- 
rable states corresponding to Ak- 
in) If we define J := N/2, 

c x := yfl «J,) 2 + (J z ) 2 )/J 2 (17) 

and p := [1 + (J x )/(Jc x )]/2 and if then Np is an integer, 
then there is also a separable state corresponding to B x . 
Similar statements hold for B y and B z . Note that this 

condition is always fulfilled, if J = and N is even. 
(Hi) There are always separable states corresponding to 
points B' k such that their distance from Bk is smaller 
than j . In the limit N — > oo for a fixed normalized angu- 
lar momentum j := J /(N/2), the difference between the 
volume of polytope of Eqs. and the volume of set of 
points corresponding to separable states decreases with N 
at least as AV/V cx A~ 2 , hence in the macroscopic limit 
the characterization is complete. 

Proof. A separable state corresponding to A. x is 

PAm :=p(\TP + )(i, + \)® N + (l-p)(\i>-)(^\f N . (18) 

Here |i/>+/_) are the single qubit states with Bloch vector 
coordinates ({a x ),(o- y ),(a z )) = (±c x , (J y ) / J, (J z ) / J). If 
M := Np is an integer, we can also define the state cor- 
responding to the point B x as 

|0 S J:=|^ M ®|V-) 0(A, - M) . (19) 

Since there is a separable state for each extreme point 
of the polytope, for any internal point a corresponding 
separable state can be obtained by mixing the states cor- 
responding to the extreme points. It is instructive to 
demonstrate this through a simple numerical experiment. 
Fig.[2]shows that for N = 10 and J = random separable 
states indeed fill the polytope. 

If M is not an integer, we can approximate B x by 
taking m := M — s as the largest integer smaller than M, 
defining 

pi := {l-e){\i, + )(^ + \r m ®{\^)(^\r {N - m) 
+ £ (|V + )(^+|)® (m+1) ® (|V-)(V-l)® (JV ~ m ~ 1) - 

(20) 

This state has the same coordinates as B x , except for the 
value of {J 2 ), where the difference is c 2 (e — e 2 ) < 1/4. The 
dependence of AV/V on N can be studied by considering 
the polytopes in the ((J 2 ), {Jy}, (J|))-space correspond- 
ing to ( Jk) — jk x yj where jk are the normalized angular 
momentum coordinates. As N increases, the distance of 
the points Ak to Bk scales as A 2 , hence the volume of 



the polytope increases as N 6 . The difference between the 
polytope and the points corresponding to separable states 
scales like the surface of the polytope, hence as A 4 . □ 

Let us analyze now our optimal spin squeezing inequal- 
ities one by one and define the corresponding facets of the 
polytope on Fig. [TJa) . Eq. (|7a[) corresponds to the facet 
A x — A y — A z . As we discussed it is valid for all quantum 
states. The symmetric states correspond to states on this 
facet and saturate Eq. (|7a|) . 

Eq. (|7b|) has already been presented in Ref. It 
corresponds to the facet B x — B y — B z . For even A, it 
is maximally violated by many-body singlets. For such 
states 

J := (0,0,0), 

K := (0,0,0). (21) 

That is, singlet states are states for which both the angu- 
lar momentum components and their variances are zero 
(28|. For large enough A there are many states of this 
type. If we mix these states, the mixture still maximally 
violates this inequality and thus it is detected as entan- 
gled. This might be the reason that this criterion can 
detect states that are very weakly entangled in the sense 
that they are separable with respect to all bipartitions. 

The violation of the criterion gives information about 
the number of spins that are unentangled with the rest 
in the following sense [13]. Let us consider a pure state 
for which the first M qubits are not entangled with other 
qubits while the rest of the qubits are entangled with 
each other 

\^) = {®lti\i>k))®\i>)M+i,...,N- (22) 

For such a state, based on the theor y o f entanglement 
detection with uncertainties, we have [291 ] 

(AJ X ) 2 + (AJ y ) 2 + (AJ Z ) 2 > M. (23) 

Let us consider now a mixed state p := J^k Pk\^k)(^k\- If 
it violates Eq. (|2"3"|) then at least one of the components 
|^'fe)( x I , fc| must have M or more spins that are entan- 
gled with other spins. If the left-hand side of Eq. (j2"3")l 
is smaller than | then the state cannot be created by 
mixing states that have one or more unentangled spins. 

Eq. ([7c|) corresponds to the facets A y — A z — B x , A x — 
A z — B y , and A x — A y — B z . All entangled symmetric 
Dicke states violate this criterion [3(|. This can be seen 
as follows. An A-qubit symmetric Dicke state with m 
excitations is defined as [3l| 

\m,N):= f^Q 2 ^P fc (|li,l 2! ...,l m ,0 m+ i,...,0 w )), 

(24) 

where {Pk} is the set of all distinct permutations of the 
spins. |1,A) is the well known A-qubit W state. For 
states of the form Eq. (|2"4"|) 

J = (0,0,m-f), 

£ = [W + !2 (^)^ + I2(A^) ;(to _^ ) 2] (25) 
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Using Eqs. (|25|) one finds that Eq. ([7c| is violated by 
all Dicke states expect for the non-entangled ones with 
m = and m = N. For even N, it is maximally violated 
by the symmetric Dicke state | tt , N). 

Finally, Eq. (|7d[) corresponds to the facets A y — B z — 
B x , A x — B z — By. and A z — B x — B y . Note that these 
inequalities detect the singlet state with Eq. ([21]) as en- 
tangled. 

Now we can ask the question, what happens if we only 
know K from Eq. ([6]) and not J from Eq. ([5]). Can we 
construct a polytope of the separable states similar to 
Observation 1? Similarly, we can consider the case that 
we know the variances [(AJ X ) 2 , (AJ y ) 2 , (AJ 2 ) 2 ], but not 
J. The following observation gives the answer. 

Observation 3. (i) Let us consider the set of points 
corresponding to separable states for even N in the 
((J x ) , (Jy) i (Jz))~ s P ace without constraining the value 
of J. This set is the polytope from Observation 1 for 
J = 0, also shown in Fig. UK®)- 

(ii) Also, the set of points corresponding to separable 
states in the [(A J x ) 2 , (A J y ) 2 , (A J z ) 2 ]-space is the same 
polytope. That is, Fig.\^a) gives also the right polytope 
if the labels of the axes are changed from (J 2 ) to (A J;) 2 . 

Proof. For the first part, it can be directly seen that 
Eqs. ([7]) are least restrictive for J = 0, for other J the 
polytope is strictly smaller. For the second part, note 
that based on Eqs. ([Jj the points corresponding to sep- 
arable states must be within the same polytope shown 
in Fig. [Tfa), even if we change the labels from (J 2 ) to 
(AJ;) 2 . It is not clear, however, that the set of separable 
states is convex in the [(A J x ) 2 , (A J y ) 2 , (A J 2 ) 2 ]-space. 
Thus, we have to show that for each separable state p 
with (J 2 ) = Si for I = x,y,z, there is a separable state 
p for which (AJ;) 2 = Si. Let us use the decomposition 
P = Y^PkPk where p k = p^ (gi pf' (g> ... <g> p)^' are prod- 
uct states. Then, such a p := ^2pkPk can be obtained by 
mixing 

Pk ■= \ {pk + JxPkJx + JyPkJy + JzPkJz)- (26) 

The state p has the same (Jf) as p. However, the value 
of (J/) 2 is zero, hence (AJ;)| = {Jf) p . □ 



III. OPTIMAL SPIN SQUEEZING 
INEQUALITIES FOR THE CORRELATION 
MATRIX 

We discuss some further features of our spin squeezing 
inequalities. One can ask what happens, if not only ( J k ) 
and (J 2 ) for k — x,y,z are known, but (Ji) and (J 2 ) in 
arbitrary directions i. We will now first show how to find 
the optimal directions x' , y 1 , z' to evaluate Observation 1. 

Knowledge of (J,) and (J 2 ) in arbitrary directions is 
equivalent to the knowledge of the vector J, the corre- 
lation matrix C and the covariance matrix 7, defined as 



Cu '■= \{JkJi + JiJk), 

Jki ■= C a - (J k )(Ji), (27) 

for k, I = x, y, z. When changing the coordinate system 
to x',y',z', vector J and the matrices C and 7 transform 
as J i-> O J, C <-> OCO T and 7 i-> 0-fO T where O is 
an orthogonal 3 x 3-matrix. Looking at the inequalities 
of Observation 1 one finds that the first two inequalities 
are invariant under a change of the coordinate system. 
Concerning Eq. (|Tc[) . we can reformulate it as 

(J 2 ) + (J 2 ) + {J 2 ) - f < (N - l)(AJ k ) 2 + (J 2 ). (28) 

Then, the left hand side is again invariant under rota- 
tions, and we find a violation of Eq. ([7c|) in some direction 
if the minimal eigenvalue of 

J:=(7V-1)7 + C (29) 

is smaller than Tr(C) — Similarly, we find a viola- 
tion of Eq. (|7d[) if the largest eigenvalue of X exceeds 
(N - l)Tr(7) - N(N - 2)/4. Thus, the orthogonal trans- 
formation that diagonalizes X delivers the optimal mea- 
surement directions x',y',z' 35]. 
Observation 4. We can rewrite our conditions Eqs. 
in a form that is independent from the choice of the co- 



ordinate system as 

Tr(C) < MM (30a) 

Tr( 7 ) > f , (30b) 

A mi „(X) > Tr(C)-f, (30c) 

A max (X) < (jV-l)Tr( 7 )-^M, (30d) 



where \ m in(A) and A max (^4) are the smallest and largest 
eigenvalues of matrix A, respectively. If Eqs. £?]) are vi- 
olated by a quantum state for any choice of coordinate 
axes x, y, and z then Eqs. \30\) are also violated. 

The preceding Observation shows how the optimal di- 
rections x, y, z can be chosen by diagonalizing the matrix 
X. However, if one diagonalizes X and does not find a 
violation of Eqs. ([50]) . this does not a priori imply that 
C, 7 and J are compatible with a separable state. The 
knowledge that for the diagonal X the off-diagonal entries 
vanish gives some additional information about the state, 
which may in principle be used as a signature for entan- 
glement. We will prove now, however, that this is not the 
case and that diagonalizing X and applying Eqs. (|30D is 
the best one can do if C, 7 and J are known. 

Note first that Eqs. (|30[) contain the following vari- 
ables: the three eigenvalues of X, Tr(C), and Tr(7). The 
latter two can be expressed with the trace of X, and J as 

Tr(C) = ±Tr(X) + ^\J\ 2 , 

Tr( 7 ) = ^Tr(X)-i|J| 2 . (31) 
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FIG. 3: The polytope of separable states corresponding to 
Eqs. (|32p for the eigenvalues of X for N — 10 and for J = 0. 
Compare with Fig. Q^a) • 



In this way, Eqs. (|30[) can be rewritten with the eigenval- 
ues of X and \J\ 2 as 



< 



N 2 (N+2) 



Tr(X) . 
Tr(X) > 4^ + \J\ 
A min (X) > ^Tr(X) 



-(N-1)\J\ S 



•^max(^) 



< 



N 
N-l 
N 



N-l 
N 



J 2 - 



Tr(X) - 



N-l 

N 



N(N-2) 



(32a) 
(32b) 
(32c) 
(32d) 



For fixed \J\ these equations describe a polytope in the 
space of the three eigenvalues of X. The polytope is 
shown in Fig. [3] The coordinates of the extreme points 
in the (Ai, A2, A3) space of the eigenvalues of X are 



N 
4 



2 N± 
k/ i 4 



El 

4 



(33) 



and 



El T \2 N 2 N 2 



(34) 



The other a/- and bk points can be obtained by trivial 
relabeling the coordinates. 

Let us now show that in the large N limit for any X and 
J fulfilling Eqs. (|52"j) there is a corresponding quantum 
state. This would mean than the conditions with X and 
J are complete and there is not another condition that 
could detect more entangled states based on knowing X 
and J. 

First, let us consider the case when J and N fulfill 
the conditions for completeness from Observation 2(ii), 
and there are quantum states corresponding to and 



bk- The states corresponding to a x and b x are pA x and 
Pb x ■= \<l>B a ){4 > B !C I, respectively, defined in Eqs. (|18H19p . 
They are the same states that correspond to the points 
A x and B x in Fig. [T] The states corresponding to the 
other extreme points can be obtained straightforwardly 
from these formulas by relabeling the coordinates. Note 
that all these states have a diagonal X matrix. Now, let 
us take a X that fulfills Eqs. (f3"2")) and diagonalize it, and 
denote it Xd after the diagonalization. Then, it is clear, 
that Xd can be obtained by "mixing" the X matrices 
corresponding to ak and ak as 



Xd 



E 



P1X1, 



(35) 



-a x ,a. 



,b x ,b*j.b z 



where pi > and ~}2 k pi — 1. Note that "mixing" X 
matrices is in general not equivalent to mixing the states, 
since X is a nonlinear function of the state. However, 
for all the states corresponding to a k and bk the vector 
J is the same and that all have diagonal X matrices. 
Therefore, the corresponding state is 



PD 



E 



pi pi- 



rn 



l=A- c ,A v .A z ,B :c .B v .B^ 



Then, if pu is the quantum state corresponding to J and 
Xd, then the quantum state corresponding to J and X 
can be obtained from pjj with coordinate rotations. Fi- 
nally, if J and N are such that no quantum state exists 
that corresponds to some of the points, then an argument 
similar to the one in Observation 2 can be applied show- 
ing that at least there is a quantum state corresponding 
to a point close to all b' k s and because of that in the 
macroscopic limit the characterization is complete even 
in this case. Thus, we can state: 

Observation 5. The criteria from Eqs. H3U\) are com- 
plete in the sense that under the conditions of Observa- 
tion 2 (ii) or for large N they detect all entangled states 
that can be detected knowing J and the correlation matrix 
C. 



IV. COMPARISON WITH OTHER SPIN 
SQUEEZING CRITERIA 

In this section we compare the optimal spin squeez- 
ing inequalities Eqs. ^\ to other spin squeezing criteria. 
First, let us consider the original spin squeezing crite- 
rion Eq. Q. This inequality is satisfied by all points 
Ak and Bk, for B z even equality holds. It is instruc- 
tive to compare the region detected by Eq. §5§ to the 
region detected by the optimal spin squeezing inequali- 
ties in the ((J|>, {J 2 ), {J 2 )) space. For a fixed J, Eq. © 
corresponds to a horizontal plane in this space, shown 
in Fig. 0Ja). Eq. ([3]) can be expressed in a way that is 
independent from the choice of the coordinate system 



Anxin(X) > |J|' 



(37) 
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Eq. ((37)) is violated if Eq. ([3]) is violated for an optimal 
choice of coordinate axes x, y, and z. 

For a state of many particles that has almost a max- 
imal spin in some direction, the standard spin squeez- 
ing inequality Eq. ([3]) is equivalent to our optimal spin 
squeezing inequality Eq. (|7c|). To see that, let us now 
rewrite Eq. (l7c|) as 



(AJ Z 



> 



1 



N 



{JD + (J 2 y) -N-l 2(N- 1)((J 2 > + (J 2 )) ■ 



(38) 



This can be transformed into 



(AJ Z ) 2 



> 



1 



N 



N-l 
(J 2 ) + 



2(iV-l)«Jj) + (J|)) 



(J 2 v) 



(J,) 2 



(39) 



Let us assume that N is large and the state has a large 
spin pointing to the x direction, that is, (Jj) ~ ^- and 
(J 2 ) « -j-. In this case (Jf) w (Jk) 2 , and the right hand 
side of Eq. (|39|) is very close to ■h. At this point one can 
recognize Eq. 

Ref. [22| presented a generalized spin squeezing in- 
equality for entanglement detection that is identical to 
Eq. (|7b[) of the optimal spin squeezing inequalities. This 
inequality has been connected to susceptibility measure- 
ments in solid state systems [HI, H(| . 

Refs. [U, [24[ presented another generalized spin 
squeezing inequality, for detecting two-qubit entangle- 
ment. According to this criterion, for states with a sep- 
arable two-qubit density matrix 



25>, 



CM N 
~ 3 



20- 



15- 



10- 




20 



<J 2 ) 40 
y 



10 



(a) 



<j2) 30 

y 

(b) 



20 



<J 2 > 30 
x 




20 ( J 2 ) 

X 



((4)+ <j?>- 



t) - 



(A-l) 2 (J m } 2 <[(^ l ) + ^ 1 ] 2 

(40) 

holds. This inequality is satisfied by all points Ak and 
.Bfc, while, when we choose (fc, i, m) — (x, y, z), for A x and 
A y even equality holds. Fig. |4jb) shows the polytope of 
the optimal spin squeezing inequality together with the 
plane corresponding to the Eq. (|40p criterion. Any state 
below the plane is detected as two-qubit entangled by 
Eq. (|40|) . Note that Eq. (|7c|) of the optimal spin squeezing 
inequalities detects all states detected by Eq. P0)l . Note, 
however, that Eq. ([40]) detects only states with two-qubit 
entanglement while Eq. (|7c)) detects entangled states that 
can have separable two-qubit density matrices. Eq. ([40]) 
can be expressed in a coordinate system independent way 



as 



A max {[^ 2 + 1 - 2Tt(C)]C - (N - 1) 2 7 } 



< [^4^)] 2 ^[Tr(C) 



Nl2 



(41) 



For states of the symmetric subspace, Eq. (1401) can be 
simplified to [H S3 



4(AJ„) 2 > i 4(J„) 2 



N 



N 2 



(42) 




20 <J 2 ) 

x 



FIG. 4: (a) Comparison of the optimal spin squeezing in- 
equalities and original spin squeezing for (k,l,m) = (x,y,z), 
N = 10 and J = (1,0,2). States detected by the latter 
are below the horizontal plane, (b) Optimal spin squeez- 
ing inequalities and the inequality Eq. (|40|l for iV = 10 and 
J = (0, 0, 0). (c) Optimal spin squeezing and criterion Eq. (|45[) 
for J= (0,0,0). 
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Violation of Eq. (|42]) for some coordinate axis z is a neces- 
sary and sufficient condition for two-qubit entanglement 
for symmetric states [33] • It can also be expressed in a 
form that is independent of the choice of coordinate axes 

m 

Kunh+TjJ^) > f- ( 43 ) 

This can be rewritten with X as 

A min (X) > (44) 

Finally, Refs. [IHjHH present a generalized spin squeez- 
ing inequality detecting entanglement close to symmetric 
Dicke states with (J z ) — 0. For separable states we have 

(J J) + (J 2) < ffiii. (45 ) 

The inequality is satisfied by all points Ak and for 
A x and A y even equality holds. Fig. [4|c) shows the poly- 
tope of the optimal spin squeezing inequality together 
with the plane corresponding to the Eq. (|45|) criterion. 
Any state corresponding to points on the right-hand side 
of the vertical plane is detected by Eq. (|45|) as entan- 
gled. Eq. (|45|) can be rewritten in a coordinate system 
independent way as 

A mi „(C) > Tr(C) - ^±ii. (46) 

V. TWO-QUBIT ENTANGLEMENT VS. 
MULTIPARTITE ENTANGLEMENT 



Interestingly, it turns our that it is also possible that 
Pav2 is separable, however, there is not an V-qubit sep- 
arable state that has p aV 2 as reduced state. In this case, 
we can conclude that the system is an entangled state 
even if p aV 2 is separable. A similar phenomenon can be 
observed in the theory of cluster states [4(|: These are 
states that are defined as eigenstates of quasi-local op- 
erators. The total state is uniquely determined by these 
quasi-local properties of the reduced states, and it can 
happen that the reduced states are separable, while the 
total state is highly entangled [4l| . 

Let us elaborate this point a little bit more. If p a v2 is 
separable and it is in the symmetric subspace then it can 
always be written in the form (2lj | 

Pav2 = ^PkPk <8> Rk (48) 

k 

with pure pu matrices. In this case, there always can be 
found an iV-qubit separable state that has this state as 
its reduced state 

PN = $>pT- ( 49 ) 

k 

Hence, we can make the following statement. 
Observation 6. // any generalized- spin squeezing 
criterion (i.e., an entanglement criterion containing 
only (Jk) and {Jk)) is violated by a symmetric state 
then the state is entangled and the reduced two-qubit 
density matrix of the state is also entangled. Note that 
this has already been known for the case of the original 
spin squeezing criterion Eq. (0) f2J/. 



Next, it is interesting to ask what kind of entanglement 
is detected by our criteria knowing that they contain only 
two-body correlation terms of the form (<r{ cr^ ) and do 
not depend on higher order correlations. In fact, all quan- 
tities in our inequalities can be evaluated based on the 
knowledge of the average two-qubit density matrix 

Pav2 := jv(^-i) E Pi i ' (47) 

where p+j is the reduced density matrix of qubits i and 
j. Do our criteria simply detect entanglement of the two- 
qubit reduced state of the density matrix? It will turn 
out that our criteria can detect entangled states with 
separable two-qubit density matrices. 

Our entanglement detection scheme is related to the 
./V-representability problem 38], i.e., to the problem of 
finding multipartite quantum states that have a given 
set of states as reduced states [3!| . When detecting en- 
tanglement based on p aV 2, we ask: Is there a separable 
N— qubit state that has p aV 2 as the average two-qubit 
reduced state. If the answer is no then we know that the 
system is in an entangled state. Clearly, if p av2 is en- 
tangled then there is not an iV-qubit separable quantum 
state that has it as a reduced state. 



On the other hand, if p aV 2 is not symmetric then it 
is possible that p aV 2 is separable but there is not an 
N— qubit symmetric separable state that has it as its re- 
duced two-qubit state. This is the reason that entangle- 
ment conditions based on J and K can detect entangled 
states with a separable two-qubit density matrix. Next, 
we will examine all our entanglement criteria from this 
point of view. 

First, we will rewrite Eqs. (|7|) as equations for the aver- 
age two-qubit density matrix p a V 2 ■ All expectation values 
in these equations are computed for p aV 2- 

S < 1, (50a) 

£ > -tt=i+ t£t E <^®!} 2 ' ( 50b ) 

i=x,y,z 

S < l + v(V m ®cr m )-(cr m ®l} 2 ^, (50c) 

£ > - 1 v l t + i^t(k®i} 2 + (^®i} 2 ) 

+ wh{v m ®or m ), (50d) 

where 

S := (a x <g> a x + a y ® a y + a z (& a z ). (51) 



9 



The first of the four optimal spin squeezing inequal- 
ities, Eq. (|7a|) , corresponds to Eq. (|50a[) . Again, this 
inequality is valid for any quantum state and cannot be 
violated. 

The second inequality, Eq. (|Tb|) corresponds to 
Eq. (|50bjl . It is known that Eq. (|7bj) can detect states 
that have separable two-qubit density matrices as have 
been shown in Ref. [271] . Such a state is, for example, one 
of the many-body singlet states, 



p s oc lim exp ( 

T->0 v 



jy 



T 



(52) 



For this state ( J") = for any m,n > 0. For such a 
state for increasing N the average two-qubit density ma- 
trix pav2 becomes arbitrarily close to the totally mixed 
state. Thus, Eq. (|50b[) is not a condition for detecting 
the entanglement of p aV 2- Moreover, note that Eq. (|7b|) 
can even detect states that are separable with respect to 
all bipartitions [26| . 

The third inequality Eq. (|7c| corresponds to Eq. (|50cjl . 
Let us consider the state 



p oc exp i 



/2 



(53) 



for N — 8 and T = 3. Direct calculation shows that 
this state is detected by Eq. ([7c]) for (k,l,m) — (x,y,z). 
Thus, again, this is not a condition for the separability 
of the two-qubit density matrix. 

The fourth condition is Eq. (|7d[) which corresponds to 
Eq. (|50dp . It detects the singlet state p s . This state has 
a separable two-qubit density matrix thus Eq. (|7dp is not 
a condition on the separability of the reduced density 
matrix. 

Let us consider now the original spin squeezing inequal- 
ity Eq. (|3|) . It is known that the violation of this inequal- 
ity implies two-qubit entanglement for symmetric states 
|42j . However, if the quantum state is not symmetric, 
Eq. ((3]) can detect states with separable two-qubit den- 
sity matrices. For example, the following state violates 
Eq. while it does not have two-qubit entanglement: 



p sq cx exp I 



2Jx Jz 

T 



(54) 



for TV = 8 and T = 0.3. 

Finally, let us consider the generalized spin squeezing 
inequality Eq. (|45| . It can be proved that any state vio- 
lating it has two-qubit entanglement. This is because it 
can be rewritten with expectation values computed for 
Pav2 as 



{a x (g> a x ) + (ay <g> cr y ) < 1. 



(55) 



Any two-qubit state violating this criterion is entangled 



VI. SPIN SYSTEMS GIVING VIOLATIONS FOR 
THE OPTIMAL SPIN SQUEEZING 
INEQUALITIES 

In the recent years, considerable effort has been made 
to create large scale entanglement in various physical sys- 
tems: In Bose-Einstein condensates of two-state bosonic 
atoms [8( , in optical lattices of cold two-state atoms real- 
izing the dynamics of an Ising spin chain [ll|, 0, [44| and 
in atomic clouds through the interaction with light and 
appropriately chosen measurements 0, [l4l [l5| . In the fu- 
ture, it is expected that experimenters will also engineer 
the various ground states of well known spin chains. En- 
tanglement detection in such systems were considered, 
for example, in Refs. [H, [H, b9, S3]. Note that there 
are methods available for measuring the variances of the 
collective spin components of atomic systems through in- 
teraction with light m mi . 

In the light of the experiments, we ask the question: 
Under what circumstances are our optimal spin squeez- 
ing inequalities useful for detecting entanglement in the 
sense that they outperform other spin squeezing entan- 
glement criteria? In this section, we will show that our 
entanglement criteria are especially useful in situations 
in which the state has a small or zero mean spin J and 
its reduced average two-qubit density matrix p aV 2 is sep- 
arable. 



A. Ground state of spin systems 

These will be, on the one hand, one-dimensional spin 
chains. On the other hand, we will consider spin systems 
corresponding to the completely connected graph. We 
will consider the following Hamiltonians: 

First let us consider the Heisenberg chain with the 
Hamiltonian 



EW +1) 



V V 



(56) 



Its ground state is a many-body singlet state. Thus, the 
optimal spin squeezing inequality Eq. (|7b[) is ideal for its 
detection. Concerning how other criteria can detect its 
ground state as entangled, we can state the following: 
Observation 7. The T = ground state of a spin sys- 
tem with a Hamiltonian without an external field can- 
not be detected by the original spin squeezing criterion 
Eq. (0). The ground state of a spin chain Hamilto- 
nian without an external field cannot be detected by the 
Korbicz-Cirac-Lewenstein criterion Eq. 
Proof. The first statement is true since criterion Eq. (|3"|) 
cannot be used for states with J = since in Eq. © one 
has to divide with the length of the collective spin com- 
ponents. The other claim can be proved noting that for 
large N the two-qubit density matrix p aV 2 of the ground 
state of spin chains without an external field is unentan- 
gled. This can be seen realizing that for the ground state 
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of an iV-qubit translationally invariant chain 
Pav2 = Tv^T (pu + Pia + Pu + ■■■ + Pin) 



(57) 



where pki is the reduced two-qubit matrix of spins k and 
I. However, for a spin chain distant sites are less and 
less correlated thus for large enough k we have pik ~ 
jl. Hence, for large enough N the reduced two-qubit 
matrix /? aV 2 is very close to the totally mixed state and 
it is separable [50]; thus, the state is not detected by the 
Korbicz-Cirac-Lewenstein criterion Eq. (pfO"|) . □ 
The Hamiltonian of the isotropic XY chain is 



H 



XY •- 



5>. 



u v u v 



(58) 



This system is similar to Eq. (|56p from the point of view 
of detecting its ground state by various entanglement cri- 
teria. That is, p av2 is unentangled for this system and the 
spin squeezing criterion Eq. © cannot detect its ground 
state. Moreover, its ground state is detected by the op- 
timal spin squeezing inequalities. While the XY chain 
is exactly solvable [5l|, the latter statement can be un- 
derstood based on simpler arguments using only quali- 
tative properties of the ground state. Let us consider a 
chain with a periodic boundary condition. For the non- 
degenerate ground state of the XY chain for even N one 
has (J|) = since J z commutes with Hxy- The near- 
est neighbor correlation is the strongest, that is for the 
ground state 



{o\ m) a[ n) ) = (-1)' 



c l,D(m,n) 



for I = x,y where D(m, n) is the distance of qubit m and 
n, and Q. m > is a monotonous decreasing function of 
to. Hence, due to translational invariance, it follows that 



N 



1 



W> = T + 7 



4 



;<7 M a („) 



V N N A 



(60) 



for l = x,y, where A N := |(ct z (1) ct, (2 } )| - | {ct z (1) cr i (3) } | . Note 
that An converges to a non-zero value for N — ► oo. Using 
these arguments, one can see that for any even N the 
ground state of the XY chain violates Eq. (|Tb[) and this 
violation is of order N in the large N limit, that is, the 
relative violation does not approach zero with increasing 
N. Hence it also follows that chains with odd N must 
also violate Eq. (|7b[) in this limit. 
The Hamiltonian 



Hr 



J T + J., 



J* 



N 

T 



=x,y 7 z m>n 



(m) (n) 
I °l 



(61) 



corresponds to a system that has a Heisenberg interaction 
between all spin pairs and has a very degenerate ground 
state. The two-qubit density matrix of its T — ther- 
mal ground state converges to the completely mixed state 
as N increases, thus for large enough N it is separable 
[27j . With respect to other qualitative statements about 
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+ 


+ 


A > 
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+ 
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TABLE I: Table showing for several spin Hamiltonians which 
entanglement condition can detect their T = thermal 
ground state in the large particle number limit. For the 
Hamiltonians see text. For the Lipkin-Meshkov-Glick model 
7 = 1 is assumed. 



entanglement detection, it is similar to the Heisenberg 
chain. 

The Hamiltonian of the Lipkin-Meshkov-Glick model 

is m 



LMG 



'N 



(J*+ 7 Jy 2 )-hJ z 



(62) 



For A > 0, 7 = 1 and h = the ground state is an 
(59) iV-qubit symmetric Dicke states with ^ excitations. For 



h all the symmetric Dicke states given in Eq. can 
be obtained as ground states of the system. These, except 
for the trivial |0,iV) = |0000...) and \N,N) = |1111...) 
states, all have entangled reduced two-qubit density ma- 
trix. Using Eq. (f2"5"l) . one can show that they are detected 
both by our optimal spin squeezing inequalities and the 
Korbicz-Cirac-Lewenstein criterion Eq. pi))) . However, 
they are not detected by the original spin-squeezing in- 
equality as can be seen by substituting Eqs. (|2"5j) into the 
original spin-squeezing inequality. For A < 0, 7 = 1 and 
h = the ground state is the same as for the Hamiltonian 
Eq. (EU). 

Finally, the summary of the results in this section is 
shown in Table |U 



B. Bound entanglement in spin chains 

Next, we study spin models in thermal equilibrium. 
We give the threshold temperatures for various spin mod- 
els for the PPT criterion [53| and for our optimal spin 
squeezing inequalities Eqs. ([7])- These temperatures are 
defined as the values, below which the spin squeezing in- 
equalities are violated of the state becomes NPT with 
respect to at least one partition. The results are given 
in Table [TTJ The systems considered are the Heisen- 
berg chain and the XY chain defined in Eqs. (|56I58[) . 
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N 


3 


4 


5 


6 


7 


8 


9 


Heisenberg 
chain 


Eq. ((7b|) 

PPT 


5.46 
4.33 


5.77 
5.47 


5.72 
4.96 


5.73 
5.40 


5.73 
5.17 


5.73 
5.38 


5.73 
5.25 


XY 
chain 


Eq. ((23 

PPT 


3.09 
2.56 


3.48 
3.46 


3.39 
3.07 


3.41 
3.34 


3.41 
3.19 


3.41 
3.32 


3.41 
3.24 


Heisenberg 
model f.c. 


Eq. (70 
PPT 


2.73 
2.16 


3.73 
2.73 


4.73 
3.17 


5.72 
3.71 


6.72 
4.17 


7.72 
4.70 


8.72 
5.17 


XY 

model f.c. 


Eq. (O 
PPT 


1.54 
1.28 


2.08 
1.82 


2.59 
2.23 


3.10 
2.74 


3.60 
3.20 


4.11 
3.71 


4.61 
4.19 


Ising chain 
B=0.5 


Eq. (JTcJ) 

PPT 


0.67 
1.08 


0.89 
1.26 


0.55 
1.17 


0.78 
1.26 


0.50 
1.21 


0.71 
1.26 


0.46 
1.22 


Ising chain 
B=l 


Eq. (JZcJ) 
PPT 


1.22 
1.49 


1.29 
1.71 


1.14 
1.61 


1.17 
1.71 


1.10 
1.65 


1.11 
1.71 


1.08 
1.67 


Ising chain 
B=2 


Eq. ([7cJ) 
PPT 


2.01 
2.15 


1.97 
2.43 


1.90 
2.30 


1.87 
2.43 


1.85 
2.36 


1.83 
2.43 


1.82 
2.38 



TABLE II: Critical temperatures for the PPT criterion and 
Eqs. (0 for Heisenberg, XY and Ising spin chains of vari- 
ous size, and for the Heisenberg any XY systems on a fully 
connected graph. For the definitions of the Hamiltonians see 
text. 
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FIG. 5: Comparison of the critical temperatures of separabil- 
ity criteria for different site numbers in the Heisenberg- (left) 
and the XY- model (right). T c of the spin squeezing inequality 
Eq. ()7b|l (+) is higher than the critical temperatures of the 
PPT (O) criterion [H or the CCNR (o) criterion [||. 



the Heisenberg system on a fully connected graph with 
the Hamiltonian Eqs. (ffTTj) . the XY system on a fully 
connected graph with the Hamiltonian -//lmg defined in 
Eq. (|62p for h = 0,7 = 1 and A = — 1, and the antifcr- 
romagnetic Ising spin chain in a transverse field defined 
as 

fl!:=j;W) + ^ ff W. (63) 

k k 

The thermal state of the system is computed as p t h °c 
exp(— -^p) with k = 1. In many cases, the temperature 
bound for the PPT criterion is lower than for our spin 
squeezing criterion. This means that there is a temper- 
ature range in which the quantum state has a positive 
partial transpose with respect to all bipartitions while it 
is still detected as entangled. Such quantum states are 
bound entangled and since all bipartitions are PPT, no 
entanglement can be distilled from them with local op- 
erations and classical communications even if arbitrary 
number of parties are allowed to join 54]. In particular, 
the results show that Eqs. ([7]) can detect fully PPT bound 
entanglement in Heisenberg and XY chains, moreover, in 
Heisenberg and XY systems on a completely connected 
graph. Note that the bound temperature for the optimal 
spin squeezing inequalities for the Heisenberg model on a 
fully connected graph is T c w N for large N |27|. On the 
other hand, our criteria do not seem to detect fully PPT 
bound entanglement in Ising spin chains. Finally, Fig. [5] 
shows the results for the Heisenberg and XY and chains, 
together with the bounds for the c omp utable cross norm 
or realignment (CCNR) criterion j5fj. The latter is of- 
ten a good indicator of bound entanglement, however, in 
these systems it does not detect bound entanglement. 



C. Bound entanglement in a nanotubular system 

Let us finally investigate a finite system showing bound 
entanglement at high temperatures. The nanotubu- 
lar system Na2Va07 is a prominent example of a low- 
dimensional quantum magnet. The compound was syn- 
thesized in 1999 by P. Millet et al. 56], who also pro- 
vided a detailed description of its structure: Every nine 
V 4+ C>5 pyramids form a ring, by sharing edges and cor- 
ners, furthermore those rings accumulate to nanotubes 
with Na atoms located in the center of and between 
them. Due to the complex structure of this system some 
years passed, until an effective model for the exchange 
interactions could be found [57t ] . The coupling terms be- 
tween the rings are considerably smaller than the inter- 
ring coupling and therefore can be neglected in a first 
approximation. Effectively the system can be described 
as a nine site anti- ferromagnetic spin-i Heisenberg ring 
showing nearest-neighbor and next-to-nearest neighbor 
interactions. The Hamiltonian can be written as 

9 . 

H := £ ^ (k) • <? (fc+1) + • a( fc+2 > , (64) 

k=l 

with periodic boundary conditions and approximately 
homogeneous parameters for the nearest-neighbor inter- 
actions Ci = 200K and C\ = 14QK for k = 2, 3, 5, 6, 8, 9, 
while C\ = in all other cases (see Fig. [6]) . The magnetic 
susceptibility of this simplified model coincides well with 
the experimental results above a temperature of about 
10K |3. 

For the given Hamiltonian, the thermal state is entan- 
gled for low temperatures and will become separable at 
a certain point when increasing the temperature. For ev- 
ery separability criterion, a critical temperature T c can 
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6 5 

FIG. 6: Schematic picture of the Na2VsC>7 system with cou- 
pling parameters Ci, C2 of the nine spin-1 Heisenberg ring 
model. 

be found. Doing so for the spin squeezing inequalities 
shows that t he c ritical temperature of the inequalities 

Eqs. (|7d]) is = 182. 8.^ while the inequality Eq. ((7b|) 

gives T<f^ = 363.6^, the other ones do not detect any 
entanglement at all. The critical temperature of Eq. (|7b[) 
has already been known from Ref. [581 ] . where the mag- 
netic susceptibility of the system has been used as an en- 
tanglement witness, which effectively results in the same 
criterion [26l |36|. Furthermore, we have computed the 
critical temperature of the Peres-Horodecki (PPT) cri- 
terion according to all bipartite splittings, resulting in a 
maximal temperature of T C PPT = 303. 9K for the splitting 
A = {1, 3, 4, 6, 7, 9} vs. B = {2, 5, 8}. So we find a transi- 



tion from free to bound entanglement at approximately 
room temperature. 



VII. CONCLUSIONS 

We presented a family of entanglement criteria that 
detect any entangled state that can be detected based 
on the first and second moments of collective angular 
momenta. We also showed that these criteria can be ex- 
tended such that they detect all entangled states that can 
be detected based on knowing the expectation values of 
the spin components and the correlation matrix. In spite 
of that these criteria do not contain multi-qubit correla- 
tion terms, they do not merely detect the entanglement 
of the two-qubit reduced state. They can even detect 
entangled states with separable two-qubit matrix. For 
further research, it would be very interesting to extend 
our results to ensembles of particles with a higher spin, 
e.g. spin-1 particles. 
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